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Abstract 

In this work we use the algebraic Bethe ansatz to derive the general 
scalar product in the six-vertex model for generic Boltzmann weights. 
We performed this calculation using only the unitarity property, the 
Yang-Baxter algebra and the Yang-Baxter equation. We have derived 
a recurrence relation for the scalar product. The solution of this re- 
lation was written in terms of the domain wall partition functions. 
By its turn, these partition functions were also obtained for generic 
Boltzmann weights, which provided us with an explicit expression for 
the general scalar product. 
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1 Introduction 



It is well-known that the quantum inverse scattering method is a powerful 
tool to solve exactly quantum integrable models as well as their classical coun- 
terparts [I1I2]- Based on this approach one can construct the Bethe states by 
the action of pseudo-particle creation operators on a pseudo-vacuum state. 
These operators are the off-diagonal matrix elements of the monodromy ma- 
trix 7^ (A). By its turn the monodromy matrix elements are the generators 
of the Yang-Baxter algebra given by the following quadratic relation, 

i?i2(A,/x)ri(A)r2(/i) = r2(^)ri(A)i?i2(A,/i). (i) 

The monodromy matrix elements act on the space of states of the quan- 
tum physical system and the i?-matrix elements play the role of the structure 
constants of the above Yang-Baxter algebra ([1]). The latter algebra is asso- 
ciative thanks to the Yang-Baxter equation [3] 

i?12(A,^)i?13(A,7)i?23(^,7) = i?23(/",7)^13(A,7)i?12(A,^). (2) 

The trace of the monodromy matrix over the auxiliary space T{X) = 
Tr^ [734(A)] is the row-to-row transfer matrix. This matrix constitutes a 
family of commuting operators [T(A), T(/i)] = 0, provided that the i?-matrix 
is invertible. This condition is granted by the unitarity relation 

i?2l(A,/x)i?i2(^,A) = /, (3) 

where / is the identity matrix. 

Taking the action of the transfer matrix over the Bethe states, one is able 
to obtain the eigenvalues of the transfer matrix and its related conserved 
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quantities. The eigenvalue expression is usually given in terms of some pa- 
rameters which should satisfy the associated Bethe ansatz equations [1]. 

After obtaining the eigenvalues, one of the most challenging problems is 
to compute the scalar products and correlation functions. This was accom- 
plished, by means of the algebraic Bethe ansatz, for models descending from 
the i?- matrix of the symmetric six- vertex model [SI [6l [7] . 

Recently, it was argued that the algebraic formulation of the Bethe states 
of the transfer matrix can be done using only the Yang-Baxter algebra ([T]), 
the Yang-Baxter equation ([2]) and the unitarity property ([3]) satisfied by the 
i?- matrix |8]. This way the on-shell properties (eigenvalues and the Bethe 
ansatz equations) as well as the off-shell Bethe vectors are obtained in terms 
of arbitrary Boltzmann weights. 

This idea have been shown to be valid also on the computation of the 
monodromy matrix elements of the arbitrary six-vertex model in the F-basis 
and the scalar product for arbitrary Boltzmann weights [S]. 

In this paper we shall obtain the scalar product as well as its associated 
domain wall partition functions by means of the algebraic Bethe ansatz for 
the completely asymmetric six- vertex model for arbitrary Boltzmann weights. 
The aim of this paper is to show that this can be accomplished using only the 
Yang-Baxter algebra, the Yang-Baxter equation and the unitarity relation. 

The outline of the article is as follows. In section [21 we define the asym- 
metric six-vertex model and list the main relations needed in this work. In 
section |3l we derive a recurrence relation for the scalar product and sketch 
its solution. In section HI we derive the partition function with domain wall 
boundary conditions which will help on writing a closed formula for the scalar 



4 



product. Our conclusions are given on sectional 



2 The six-vertex model 



(4) 



We define the i?-matrix of tlie asymmetric six- vertex model as [31 [ID] , 

^ a+(A,^) ^ 

6+(A,/i) c+(A,/i) 

c_(A,/i) 6_(A,/i) 

\ a_(A,/i) y 

This i?-matrix should satisfy the unitarity property which results on the 
following set of relations for its matrix elements 



c_(A,/i)6„(/i, A) - 


f 6_(A,/i)c+(/i, A) 


= 0, 


6+(A, fi)c_{ij,, A) + c+(A, A) 


= 0, 




a+(A,/i)a+(/i, A) 


= 1, 




a_(A, fi)a^{fi, A) 


= 1, 


c_(A,/i)c_(/i, A) - 


f 6_(A,/i)b+(/i, A) 


= 1, 


c+(A,/i)c+(/i, A) - 


f 6+(A,/i)6_(/i, A) 


= 1. 



(5) 
(6) 
(7) 
(8) 
(9) 
(10) 



Moreover the above i?- matrix satisfies the Yang-Baxter equation ([2]). This 
provides us with an additional set of functional relations among the Boltz- 
mann weights given by, 



o_c„o , + c_a , c_ 



6+a' c': + c+c'_b'' 



Cl-|_C_Cl_|_ , 



(11) 
(12) 

(13) 







(14) 


h+d_h"_ + c-a'_d'_ 


= a_c_a_, 


(15) 


b.a'_c'L + c+c'^b'L 


= a_6'_c':, 


(16) 




= a+f)'_,_c'|, 


(17) 


0-C_^_b_^_ + c+a_,_c_,_ 


/ // 

— C_j_ (X_j_ , 


(18) 


c^a',?)'; + 6+c'_c'| 




(19) 




= a^b'^cl, 


(20) 


c+a_b^ + o+c_,^c_ 


1,1 II 
= c+o_,_a_, 


(21) 




= a-c'_^_a'!_, 


(22) 



where a± = a±(A,//), = a±(A,7) and a'^ = a±(/^,7) and likewise for the 
other functions. 

The monodromy matrix 7^(A, {i^k}) = Ral{^, ^l) ■ ■ ■ -R^i(A, i^i) satisfies 
the fundamental relation ([T]) thanks to the Yang-Baxter equation ([2]). For 
the six- vertex model, the monodromy matrix can be written as a 2 x 2 matrix 
with operator valued entries, 



r(A, {/.,}) 



(23) 



C{\,{u,}) DiX,{jyk}) 

These operators have to obey sixteen commutation rules which follow from 
([1]). We only list the commutation rules used in this work. 



a+(A,/u) 



(24) 



6„(/i, A) 



6_(/i,A)' 



mSif^) = ?4M^(/^)^W - ^#4^W^(/^)> (26) 



6„(A,/i) 



6„(A,/i)- 



c{X)C{^,) = ^4^^(^)C(A), (27) 

a+(/i, A) 

Additionally, the structure of the i?-matrix (jlj) implies that both ferro- 

magnetic states |ft) = (S),=i ^^(^ l-ll) = 01=1 innate 

eigenstates of the transfer matrix T{X,{uk}) = 74(A, {z/fc}) + D{X, {uk}). 
Therefore, they play the role of pseudo-vacuum states. 

One can see immediately that the action of the monodromy matrix over 
the pseudo-vacuum state, e.g. Ifl'), results in a triangular matrix, 

/ 



(31) 



aiX,{i^,})m # 
y d{X,{uk})\r, 

where a(A, {i^k}) and d{X, {I'k}) correspond to some fixed representation and 
the symbol # stands for a non-null state. Therefore, one sees that the 
B{X,{uk}) and C(A,{i^fc}) operators are the creation and annihilation op- 
erators over the state Ift) |4], which means 

^(A, M) m = a(A, {uk}) m B{X, {//,.}) 1^) = #, (32) 
C(A, {lyk}) m = D{X, {u,}) \t) = d{X, {ly,}) m , (33) 

or alternatively, 

(^1 A(A, {z/fc}) = (^1 a(A, {uk}) m BiX, {u,}) = 0, (34) 
{t\ C{X, {u,}) = # (^1 D(A, {u,}) = {t\ d{X, {//,}). (35) 
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The action of the creation operators over the pseudo- vacuum is the alge- 
braic version of the famous Bethe ansatz and is given by 



(36) 



In order to obtain the eigenvalues of the transfer matrix one has to com- 
mute the operators A{X, {I'k}) and D{\, {z^fe}) with B{\, {^k})- In doing so 
we need to use the relations fl20|24|25l) for the operator A{X, {lyk}) and the re- 
lations fll2|24|26p for the operator D{X, {t'fc})!!], which results (see appendix 
A) 



A{x, M) n = n rjY^ n ^(^- m)a{x, {^.d 



i=l 



;^&-(A.A)t^ 



g n ^^^(^. {-}) n M), (37) 



M 



M 



i = l 

M 



D{x, n Bik, = n ?4rTT n ^(^- M)Dix, 



i=l 



^■^fe-(A,A, 



M 



i=l 
M 



where a^^(Aj,Aj) = a+(Aj, Aj)6'>(Aj, Aj) and 



^>(Aj, Aj 



g- (Ai,Aj) 
a+(Ai,Aj) 

1, 



(39) 



Therefore, we obtain 



M M 

T(A) |M/m) = Aa/(A, [vk]) I^A/) + 5^ r,(A)i?(A, {uk}) J] ^(A., {j^J) |^) ,(40) 



!=1 
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where the eigenvalues Aa/(A, {I'k}) are given by 



M , . M 



A.,(A. M)^ a(A, n ^ H- ^(A. {..}) R («) 
and 

Finally, we have to impose the coefficient rj(A) to vanish for arbitrary A 
in order to cancel the unwanted terms. Using the unitarity relation ([5]), one 
sees that the parameters Xj have to satisfy the Bethe ansatz equations 

«(\-^{^fc}) ^ TT Q?(Aj, A^) b-(A„Aj) . . 

-L|af(A.,A,)&-(A„A.)- ^ ^ 

Similarly the dual Bethe state can be written as follow, 

(vI/mI = (^1 C(Ai, {uk}) ■ ■ ■ C(Am, {uk}). (44) 

However, in order to obtain the transfer matrix eigenvalues one should look 
for commutation rules between the operators A{X, {t'/t}) and D{X, {I'k}) with 
C(A, {z/fc}). Here, one has to use the relations (I16f27ll28p for the operator 
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A{X, {i^k}) and the relations fll4|27|29p for the operator D{X, {j^k}), such that 

n {^^k})A{\ = n ^7x^^(^' {^4) n ^(^- {^^^^ 

i=l i=l V *' / j=i 

- E ^ f[ Ig^^t^. n ^^(A.. {^.}). (45) 
n c(A., {/..})D(A, {z.,}) = n ^ttat^^^' ^'''^^ n ^(^- {^^^i) 

i=l i=l ^ ' i=l 

3 The scalar product 

The general scalar product is defined as 

SAimfLmZi) = (^|C(/ii)---C(M5(AM)---i?(Ai)|^) (47) 

where we have dropped the dependence on i/fc. In the particular case that 
both set of parameters A^ and fij satisfy the Bethe ansatz equations fj^3|) . the 
product p7|) becomes the norm of Bethe eigenstates[5]. 

To compute the scalar product (jH]), we must calculate the action of the 
operator B{X) over the state (jUj) (or alternatively, the action of C{fi) over 
the state (!36|) ). To this purpose we have to use repeatedly the relation (!30|l . 
which results 

(^1 C(Azi) ■ ■ ■ CMB{Xm) = n . f x"" (^1 5(AM)C(/ii) ■ ■ ■ CM 

+ /I'^X) n^(/^-) (-4(AM)/^(/i,) - A(/x,)D(Am)) n C'(/x.),(48) 
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where the first term vanishes thanks to (IMI) and we set Aj = fi2M+i-i, i = 
1, ■ ■ ■ , M for convenience. 

After that, one has to compute the action of A{X)D{fij) — A{fij)D{X) 
over the C{fii), where one should use the algebraic relations fl45H46p along 
the same lines as [21 15], such that 

A/+1 

(^1 C7(/ii) ■ ■ ■ C{fiM)B{fiM+i) = Yl ^M^if^i) X (49) 

j,fc=i 

i^k i^k 



where we have also used the functional relation f[T6|) and the unitarity rela- 
tions ([3l9]) for simplifications. 

So we multiply the previous expression by additional operators, 

SAimi^M+imfii) = m ci^i) ■ ■ ■ c{^m)b{^m+i) n-iU^ ^(/^.) i^)- tws 

results in a recurrence relation for the scalar product given by. 



SMmi=M^imfii) = E 



a{fik)d{iJ,j)c_{fij, hm+i)c+{hm+i, IJ'k)a'+ ifJ'j, l^k) 

M+l (0)/ ^ (^),' N 

-7—7 \~n ^^M-l{{^J'\i=M+2\\^^\i^l)■ \P^) 



X 



One can iterate the above recursion relation ( l50ll and represent the re- 
sulting expression as follows, 

{A} = {A+}U{A-} i=l t=l J IP J 

{f} = {M+}U{M-} 

I A~t~ I + 1 M I — 
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where we sum with respect to all partitions of the set {A} U {/i} into two 
disjoint sets {A'^} U and {A~} U of M of elements. If the number 
of elements in the sets {A"^} and {i^t'^} is n+ = \\~^\ = |^^|, then we have 
?7,_ = |A~| = |//_| = M — elements in the sets {A~} and 

Notice that any coefficient Km is determined only by terms arising from 
algebraic relations among monodromy matrix elements. Therefore, Km is 
independent of the representation, i.e. independent of the choice of the func- 
tions a(A, {i^k}) and d{X, {i^k}) 

In view of that one can fix the coefficients Km using any special repre- 
sentation. Let's consider the case where L = M and 

M M 

aM(A) = JJa+(A,z/i), (ImW = Y[b4X,iy,), (52) 

i=l 1=1 

where we have fixed the inhomogeneities as follows. 



(53) 



A+, i = !,■■■ ,n+ 

Using the fact that the unitarity relation ([5]) implies that 6_(A, A) = 0, 
one sees that (^^/(A) = for any A G {A+} U and ^^/(A) 7^ for any 

A G {A^}U{/i'^}. Consequently, there is only one non- vanishing term in the 
sum flSTl) . which allow us to write 

l[aM{Xt)dM{^it)l[aM{^^i)dM{X-)KM{{X^}, {A-}|{/i+}, {/i"}) 

i=l i=l 

= {UmI C(/ii) ■ ■ • C(//m)S(Am) ■ ■ ■ 5(Ai) IHm) , (54) 
where Iff'M) is a state with M spins up. 
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The product of 5-operators overturns all M spins resulting in the state 
I-IIm). Therefore, we have 



(^m| C(Aii) ■ ■ ■ C{^im)B{Xm) ■ ■ ■ 5(Ai) IHm) = 
= zS\{f^h {A+} U ({A}; {A+} U (55) 

where the functions Z^'^^ are defined by 

4?({A};{A+}U = (^m|S(Am)---S(Ai)|^m), (56) 

^l?(M;{A+}u {//_}) = (^m|C7(/xi)---C^(MIW- (57) 

Finally, the coefficient Km can be written in terms of the above functions as 
follows 

{A+} {A+} U {^-})4f ({A}; {A+} U {^_}) 



K 



M 



(58) 



{A } {/^ }y nr=i«Af(A+)rfM(^^+)nr=iaAK^i )c^A/(Ai ) 

The functions f l56H57|) are usually called domain wall partition function 
[S] and play an important role on the calculation of scalar products and 
correlation function. We shall compute these partition functions on the next 
section in order to write a closed formula for the scalar product (!5T]) . 



4 The domain wall partition function 

In this section we derive a recurrence relation for the partition function for 
the asymmetric six-vertex model with domain wall boundary conditions for 
arbitrary Boltzmann weights. This way we will proceed along the same 
lines of [TT] and define a couple of auxiliary one-point boundary correlation 
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functions as follows, 



G'Sn = {^m\ B{Xm) ■ ■ ■ B{Xn+i)pYB{Xn)B{Xm-i) ■ ■ ■ B{Xi) , (59) 
Hm!n = B{Xm) ■ ■ ■ B{XN+i)PiB{XN)ptB{XN^i) ■ ■ ■ B{Xi) lUi) ,(60) 

where pf = ^{1 ±(tI) are the local spin up and down projectors and G^^j\,j = 
1. The above boundary correlations are related by 

N 

MB) _ tt{B) :MB) ff^r,^ 
^M,N — ^M,N ^ '^M,N-1- y^^) 

One can use the above relations in order to write the domain wall partition 
function Z^j in terms of the boundary correlations 

4f ({A}; {^}) = Yl (^^^\ ^(^^^) ■ ■ ■ B{X,^,)p^B{X,)ptB{X,^,) ■ . ■ S(Ai) .(63) 
i=i 

To derive the recurrence relation for the partition function we use the 
two-site model decomposition. More specifically, one has to decompose the 
monodromy matrix into two parts and introduce two monodromy matrices 

rA{X) = Tf{X)Ti'\X)={ ^^^^ ^^^^ I, (64) 

C{X) D{X) 



such that 



r« = i?^i(A,.o=( ^^^'^ ""'^'^ I, (65) 

Ci(A) Di(A) 

ro^ / A2(X) B2(X) , 

rj'^ = Ram{X, um) ■■■Ra2{X,U2)= \ \, (66) 



C2(A) D^iX) 
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which act on the states Ifl'j), i = 1,2 given by 

ii)= I ^ j , 1^2) = (8);i2 1^ ^ I . (67) 

In particular, the required monodromy matrix elements are readily ob- 
tained 

B{X) = A2(A)fii(A) + S2(A)Di(A), (68) 
Si(A) = c+(A,i^i)ar, (69) 
DiiX) = b.{X,i^,)pt + a.{\,ui)pi, (70) 

and the action of the operator A2{X) on the state |ff2) is given by 

M 

^2(A)|^2) = n«+(^'^'')l^2). (71) 

i=2 

Using fl68ti70p . we can rewrite (ES]) in terms of the operators A2{X) and 
i?2(A) as follows, 

M j-1 M 

ziS\{X}; {u}) = Y.c^{X,,u,)l[b_{\,u,) H a.{\,u,) (72) 

j=l i=l i=j+l 

X (^2| 52(Am) ■ ■ ■52(A,+i)A2(A,)S2(A,_i) ■ ■ ■52(Ai) 1^2) . 

At this point one should use the algebraic relations (124113 7p [TT] in order to 
obtain the recurrence relation for the domain wall partition function, 

M M M 

j=l 1=1 i=l 
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Iterating the above relation M — 1 times, one gets the solution as a sum 
over all permutations P of {X}fii 



p 1=1 I,J=1 ^ 

(C) 

We still have to determine the partition function Z^^ . The calculation 
goes in a similar way as above. Therefore, we just present the final result for 
the recurrence relation, 

M M M 

M {9), X 

whose solution obtained by iteration is given as a sum over all permutations 
P of {^^]fU 

M M (9)/ _ _ \ 

5 Explicit expression 

Using the special choice for the inhomogeneities (153|) and the recurrence 
relations fl73ll75p . we can write a simplified expression for the coefficient Km 
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1111 h {,,+ n--) h (\- \ + -) 11 



as follows 

j<k 

^ zS>{{^i^}■,{x^})zl\{x-}■,{^.-}) 

Finally, if we substitute the expression f l7i|76|77p in fl5T]) . we obtain an 
explicit expression for the scalar product for arbitrary Boltzmann weights, 

c (f\\\f w - nr=i aiK)d{fj't) nr=i q(^»~)^(a^~) 

l^+l + l)'-|="+ 



X 



-pr -j^ Q+(/^^/^fc)a-(Afc,A+) -pj. a-(A+,A+) 

M Ik ^.-) ^-(^.-' a;) 11 a+(A;, A+) ^ 

j<fc 



X 



6 Conclusion 

In this paper we obtained an explicit expression for the scalar product and 
the domain wall partition function for the arbitrary six-vertex model. 

We have reviewed the algebraic Bethe for the six-vertex model under 
the new perspective of working with arbitrary Boltzmann weights. Using 
the main ingredients of the algebraic Bethe ansatz, we have been able to 
derive a recurrence relation for the scalar product. We managed to write 
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the solution of this recurrence relation in terms of the domain wall partition 
functions. Then, we obtained a closed formula for the required partition 
functions. Finally, we have obtained an explicit expression for the general 
scalar product in terms of the arbitrary Boltzmann weights. We have done 
all that only using the Yang-Baxter algebra, the unitarity relation and the 
Yang-Baxter relations. 
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Appendix A: two and three-particle state 

In this appendix we work out in detail the commutation of the operator A(X) 
with two and three 5-operators in order to exemplify the use of Yang-Baxter 
relations flim22p . In fact, to pass the 74-operator over any number of B's one 
need only to use the Yang-Baxter relation fl20p . 
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A.l M = 2 



Let's start with the two S-operators case and use the commutation rule ([25 
twice, which results 

a+(A2, A) a+(Ai, A) 



A(A)S(A2)S(Ai) = 

c+(A2, A) a+(Ai, A2) 
b_(A2,A) 6„(Ai,A2) 
- a+(A2, A) c+(Ai, A) 
.6_(A2,A) 6-(Ai,A) 



6_(A2,A) 6_(Ai,A) 
B{X)B{X^)A{X2) 



S(A2)S(Ai)A(A) 



(A.l) 



5(A2)5(A)-^^^tH4i?(A)S(A2) A{X,). 



b_(A2,A) b_(Ai,A2)^ 

We can further manipulate the third term of the above expression using 
the commutation rule (l24l) and the unitarity relation ([5]) as indicated. This 
provides us with the following expression 
a+(A2, A) a+(Ai, A) 



A{X)B{X,)B{X^, 

6_(A2, A) 6-(Ai, A) 

c+(A2, A) a+(Ai, A2^ 



B{X2)B{X^)A{X) 



6_(A2,A) 6_(Ai,A2) 



B{X)B{Xi)A{X2) 



eg. 



b_(A2, Ai)a„(A2, A)c+(Ai, A) + c_(A2, Ai)c+(A2, A)6„(Ai, A) 



(A.2) 



S(A)fi(A2)A(Ai 



6_(A2,Ai)&_(A2,A)fc_(Ai,A) 
Finally, one can see that the resulting expression for the third term coincides 
with the left hand side of the Yang-Baxter relation fl2Ul) . After substituting 
([20D in (ICT) . we obtain 

a+(A2, A) a+(Ai, A) 



A(A)S(A2)5(Ai) 



b_(A2,A) 6_(Ai,A) 
c+(A2, A) a+(Ai, A2) 

b4X2,X) fc-(Ai,A2) 
Ch-(Ai, A) a_(A2, Ai) 
{>_(Ai,A) 6_(A2,Ai) 



5(A2)5(Ai)A(A) 
i?(A)i?(Ai)A(A2) 
B{X)B{X2)A{X^) 



(A.3) 
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A.2 M = 3 



Now we turn to the case where we have three S-operators. We use the 
previous result (1A.4I) and the commutation rule ( 12^ . such that 

a+(A3, A) a+(A2, A) a+(Ai, A) 



c+(A3, A) a+(A2, A3) a+(Ai, A3) 



6„(A3,A) 6-(A2,A3) 6-(Ai,A3) 
c+(A2, A) a_(A3, A2) a+{Xi, A2) 



+ 



6„(A2,A) 6_(A3,A2) b_(Ai,A2) 

■ _ c+(Ai, A) a+(A3, A) a+(A2, A) 
[ b4Ai,A)fo_(A3,A)6_(A2,A) 

c+(A3, A) c+(Ai, A3) a+(A2, A3) 



5(A)fi(A2)S(Ai)A(A3) 
BiX)BiX,)BiX,)A{X2) 



(A.4) 



eg. 



S(A3)5(A2)5(A) 



6_(A3,A) 6_(Ai,A3) 5_(A2,A3 
c+(A2,A) c+(Ai,A2) a_(A3,A2) 



B{X)B{X2)B{Xs) 

eg.|2l 

i?(A)5(A3)i?(A2)lA(A0. 



6_(A2,A) MAi^&-(A3,A2) 

eg.|5l 

Again one should use the commutation rule (l24l) and unitarity relation ([5]) 
as indicated. So the last term of ( ]A.4|) can be rewritten as follows, 

c+(Ai, A) a-(A3, A) a-(A2, A) _ c+(A3, A) c+(Ai, A3) a-(A2, A3) 
6_(Ai, A) 6_(A3, A) 6_(A2, A) 6„(A3, A) 6_(Ai, A3) 6-(A2, A3) 

At this point, we have to simplify the term inside the square bracket (Ji). 
This can be done by joining the first and the last term on ( 1A.5I) and using 
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the equation (1201) a couple of times, such that 

[6_(A2, Ai)a_(A2, A)c+(Ai, A)]a_(A3, A)6_(A3, A2) 



h 



+ 



6_(Ai, A)6_(A3, A)6_(A2, A)6_(A2, Ai)6_(A3, A2) 

eg. 1201 



a_(A3, A2)6-(A3, A)c+(A2, A)] c_(A2, Ai)6_(Ai, A) 



6_(Ai, A)6_(A3, A)6_(A2, A)6_(A2, Ai)6_(A3, A2) 
c+(A3, A) c+(Ai, A3) a-(A2, A3) 
6_(A3,A) 6_(Ai,A3) 6_(A2,A3)' 



(A.6) 



and 



eg 



+ 



^ a„(A3, A) [6_(A2, Ai)a_(A2, A)c+(Ai, A) + c_(A2, Ai)c+(A2, A)b_(Ai, A)] 

6-(A3,A) 6_(A2,Ai)6_(A2,A)6_(Ai,A) 
c_(A3, A2) c+(A3, A) c_(A2, Ai) c+(A3, A) c+(Ai, A3) a_(A2, A3) 



MAg^ 6„(A3, A) 6_(A2, Ai) 6„(A3, A) 6_(Ai, A3) 64A2, A3) 

eg.lSl 

Then we again have to use the equation ( l20l) twice, 
c+(Ai, A) a_(A3, A) a-(A2, Ai) 



(A.7) 



6_(Ai,A) b_(A3,A) 6_(A2,Ai) 



(A.8) 



eg. 



c+(A3, A) [c„(A2, Ai)c+(A2, A3)fe_(Ai, A3) + ^-(Aa, Ai)a_(A2, A3)c+(Ai, A3)] 
fe-(A3,A) 6_(A2,Ai)b_(A2,A3)6_(Ai,A3) 
a_(A2, Ai) r c+(Ai, A) a-(A3, A) _ c+(Ai, A3) c+(A3, A) 
{>_(A2,Ai) [6_(Ai,A)fe_(A3,A) 6_(Ai, A3) 6_(A3, A) 

eg.ini 



(A.9) 



and 



eg. 



«-(A2, Ai) [6_(A3, Ai)a_(A3, A)c+(Ai, A) + c_(A3, Ai)c+(A3, A)b_(Ai, A)] 



&™(A2,Ai 



6_(A3,Ai)6_(A3,A)6_(Ai,A) 



a-,(A2, Ai) a-(A3, AQ c+(Ai, A) 
6_(A2,Ai) 6_(A3,Ai) b4Ai,A)' 



(A.IO) 
(A.ll) 
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Finally, we can substitute the final result for the term Ii (lA.lip in Eq. (lA.4p . 

B{X)B{X2)B{X,)A{X,) 
B{X)B{Xs)B{X,)A{X2) (A.12) 
B{X)B{X3)B{X2)A{Xi). 



c+iX-s, 


A) a+(A2, A3) a+(Ai, 


,A3) 




A) 6-(A2,A3) 6_(Ai, 


A3) 


C+(A2, 


A) a_(A3, A2) a+{Xi, 


,A2) 




A) 6_(A3,A2) 6_(Ai, 


A2) 


c+(Ai, 


A) a_(A3, Ai) a_(A2, 


,Ai) 



6„(Ai,A) 6_(A3,Ai) 6-(A2,Ai) 

This expression coincides with the formula ( 1371) for the case M = 3. Note 
that we have only used the commutation rules fl24ll25p . the unitarity relation 
([5]) and the Yang-Baxter relation (120|) . It is remarkable that the same applies 
to any M-values, where we again would need only the mentioned relations. 

Similarly, the expression (l38l) can be obtained using the commutation 
rules fl24|26p . the unitarity relation ([5]) and the Yang-Baxter relation f|T2l) . 

Alternatively, one could have obtained the formula fl37j) in a shorter way. 
Note that the first term in (1371) is obtained straightforwardly using the first 
term of the algebraic relation (125!) M-times. On the other hand, the second 
term of (1371) can be very complicated, as we have seen above. The coefficient 
of the term not containing the first i?-operator, B{Xm), is very simple though. 
This coefficient is obtained using the second term of ( l25l) in order to exchange 
the A{X) and B{Xm) and in all other steps one has to use only the first term 
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of fl25p which results in the following expression 

M M M 

i=l i=l V / 

+ sum of other terms not containing B{Xj), j = 1, ■ ■ ■ , M — 1. 

However, one could have done similar analysis for any B{Xj). In doing 
so, one has to use the relation (12^ in order to move B{Xj) to the leftmost 
position, 

M 

B{\m) ■ ■ ■ B{Xj+,)B{Xj) ■ ■ ■ S(Ai) = B{X,) n B{X,)9^{X„ A,) (A.14) 
» V ' 



i=i 



eo. l|24|l i/i 

and then repeat the above described procedure. This would result in the 
missing terms in (]A.13I) and finally we should obtain the expression ( 1371) . 



References 

[1] L.A. Takhtajan and L.D. Faddeev, Russian Math. Surveys, 34 (1979) 
11 

[2] V.E. Korepin, G. Izergin and N.M. Bogoliubov, Quantum Inverse Scat- 
tering Method and Correlation Functions, Cambridge University Press, 
1993 

[3] R.J. Baxter, "Exactly Solved Models in Statistical Mechanics", Aca- 
demic Press, New York, 1982. 

23 



[4] L.A. Takhtajan, "Introduction to algebraic Bethe ansatz", Lecture 
Notes in Physics, Vol. 242, 1985, Springer- Verlag, p. 158; L.D. Fad- 
deev, "How algebraic Bethe ansatz works for integrable models", Les 
Houches School on Quantum Symmetries, 1995, North-Holland, p. 149 

[5] V.E. Korepin, Comm. Math. Phys., 86 (1982) 391 

[6] A.G. Izergin and V.E. Korepin, Commun. Math. Phys., 94 (1984) 67; 
Commun. Math. Phys., 99 (1985) 271 

[7] N. Kitanine, J.M. Maillet and V. Terras, Nucl. Phys. B, 554 (1999) 
647; Nucl. Phys. B, 567 (2000) 554; N. Kitanine, J.M. Maillet, N.A. 
Slavnov and V. Terras, Int. J. Mod. Phys. A, 19 (2004) 248 

[8] C.S. Melo and M.J. Martins, Nucl. Phys. B, 806 (2009) 567 

[9] M.J. Martins and M. Zuparic, Nucl. Phys. B, 851 (2011) 565 

[10] R.J. Baxter, J. Stat. Phys., 108 (2002) 1; Studies in Appl. Math. 
(M.I.T.) 50 (1971) 51 

[11] N.M. Bogohubov, A.G. Pronko and M.B. Zvonarev, J. Phys. A: Math. 
Gen. 35 (2002) 5525 



24 



